Two-species fermion mixtures with population imbalance 



O 

o 

(N 

< 
m 

o 
o 

u 

Oh' 



I 

O 

o 



> 

oo 

o 
o 

-(— > 



o 



X 



M. Iskin and C. A. R. Sa de Melo 
School of Physics, Georgia Institute of Technology, Atlanta, Georgia 30332, 

(Dated: February 6, 2008) 



USA 



We analyze the phase diagram of uniform superfluidity for two-species fermion mixtures from the 
Bardeen-Cooper-SchriefFer (BCS) to Bose-Einstein condensation (BEC) limit as a function of the 
scattering parameter and population imbalance. We find at zero temperature that the phase diagram 
of population imbalance versus scattering parameter is asymmetric for unequal masses, having a 
larger stability region for uniform superfluidity when the lighter fermions are in excess. In addition, 
we find topological quantum phase transitions associated with the disappearance or appearance of 
momentum space regions of zero quasiparticle energies. Lastly, near the critical temperature, we 
derive the Ginzburg-Landau equation, and show that it describes a dilute mixture of composite 
bosons and unpaired fermions in the BEC limit. 

PACS numbers: 03.75.Ss, 03.75.Hh, 05.30.Fk 



Major experimental breakthroughs have been made re- 
cently involving one-species trapped fermions (^Li) in two 
hyperfine states with different populations. The super- 
fluid to normal phase transition and the vortex state [l| , 
as well as phase separation between paired and unpaired 
fermions [2| were identified as a function of population 
imbalance and scattering parameter. These studies are 
important extensions of the so-called Bardeen-Cooper- 
Schrieffer (BCS) to Bose-Einstein condensation (BEC) 
evolution for equal populations, which were studied via 
the use of Feshbach resonances |M lj> IE IE IZ') S,J • The prob- 
lem of fermion superfluidity with population imbalance 
has been revisited recently in several theoretical works 
in continuum Ji ES O iH |ll S IM El IH 111 El 
and trapped Hll lil li3,131illii| atoms. These 
one-species experiments are ideal candidates for the ob- 
servation of uniform and non-uniform superfluid phases, 
which may be present not only in atomic, but also in 
nuclear (pairing in nuclei), astrophysics (neutron stars), 
and condensed matter (superconductors) systems. 

Arguably one of the next frontiers of exploration in 
ultracold Fermi systems is the search for superfluidity 
in two-species fermion mixtures (e.g. ^Li and '"^K) with 
and without population imbalance. While earlier works 
on two-species fermion mixtures were limited to the BCS 
limit [si llOJ I , in this manuscript we study the evolution of 
superfluidity from the BCS to the BEC limit as a function 
of the scattering parameter and population imbalance. 

Our main results are as follows. At zero temperature, 
we construct the phase diagram for equal and unequal 
masses of paired fermions as a function of scattering pa- 
rameter and population imbalance P — {N-^ — Ni)/{N'^ + 
Ni) as shown in Fig. ^ The phase diagram is asym- 
metric for unequal masses, having a larger stability re- 
gion for uniform superfluidity when the population iV| of 
lighter fermions is larger than the population N^ of heav- 
ier fermions. In addition, we find topological quantum 
phase transitions in the phase diagram associated with 
the disappearance or appearance of momentum space re- 
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FIG. 1: Phase diagram of P = {N-f - Ni)/{N^ + iVJ versus 
l/{kF,+aF) for a) equal {m-f = m^) and b) unequal masses 
(m| = O.I5m|). We show normal (N), non-uniform (NU) or 
uniform (U) superfluid phases. The dotted and P — lines 
separate topologically distinct regions. In b) the U phase also 
occurs for P < when l/{kF,+aF) > 4.8 (not shown). 



gions of zero quasiparticle energies. Lastly, near the criti- 
cal temperature, we derive the time-dependent Ginzburg- 
Landau (TDGL) equation, and show that it describes a 
dilute mixture of bosons (tightly bound fermions) and 
excess (unpaired) fermions in the BEC limit. 

To describe a dilute two-species Fermi gas in three di- 
mensions, we start from the Hamiltonian {h = 1) 



H ^ 






^^Wk 



k,^"k,^ 



E nK^')bU 



'k',q, 



(1) 



k,k',q 



where the pseudo-spin a labels the hyperfine states rep- 
resented by the creation operator aj^^, and 6J = 

"k+q/2.T"-k+q/2,r ^^^^^ Ck,o- = Ck.o- " ^J■<T, whcrc ek,o- = 
k"^ I {2m„) is the energy and /ig. is the chemical potential 
of the fermions. Notice that, we allow for the fermions to 
have different masses m„ and different populations con- 
trolled by independent chemical potentials ^a- The at- 
tractive fermion-fermion interaction l^(k, k') can be writ- 
ten in a separable form as T^(k, k') = — f/FJ^Fk' where 
5 > 0, and Fk = 1 for the s-wave symmetry. 



The gaussian effective action for H is |27| 5'gauss = 
So + (/?/2)E,At(9)F-i(g)A(g), where q = {q,Vi) with 
bosonic Matsubara frequency Vi = 2lTr/(3. Here, /3 = 
1/T, A^q) is the order parameter fluctuation field, and 
the matrix F~^{q) is the inverse fluctuation propagator. 
The saddle point action is 



5o=/3^ 



J2 {/3(^k,+ - i?k.+) 

k 

ln[np(-£;kj)] +lnK(-i?k,i)]},(2) 



where _Bk,cr = (Ck,+ + l^kP)^/^ + s^Ck,- is the quasipar- 
ticle energy when S| = 1 and the negative of the quasi- 
hole energy when S| = —1, and Ek^± = (i?k,T i ^k,i)/2- 
Here, Ak = AoFk is the order parameter, ni^(ii^k,cr) is 
the Fermi distribution and ^k,± — [S.'kA i ^k,j.)/2 — 
k'^/{2m±) — fj.±, where m± ~ 2m,fTO|/(77i| ± mi) and 
/i± = (/i| ± fJ.i)/2. Notice that m+ is twice the reduced 
mass of the t and | fermions, and that the equal mass 
case corresponds to \m-\ —^ oo. The fluctuation term in 
the action leads to a correction to the thermodynamic 
potential, which can be written as figauss = f^o + ^fiuct 
with Qo = So/P and n^^.t = (1//3) ^^lndet[F-i(g)//3]. 
The saddle point condition SSo/SA^ = leads to an 
equation for the order parameter 



9 ^ 



2E^ 



k,+ 



"-^kH 



(3) 



where Ak,± = (-^k,! i Ak_|)/2 with Ak,j = 
tanh(/3i5k,(T/2). As usual, we eliminate g in favor 
of the scattering length ap via the relation 1/g — 

-m+y/(47raF) + Ek |rk| V(2ek,+), where ek,± = (ek,T ± 
ek,x)/2. The order parameter equation has to be solved 
self-consistently with number equations N„ = —dU,/d^a 
which have two contributions N„ = N, 



0,0-" 



-A^fluct.O-- ^^0,(7 



is the saddle point number equation given by 



No,.^ 






1 - s„Xi,__ 



2-Ek + 



A'k 



(4) 



and A^fluct,<T = —dflauct/d^ia- is the fluctua- 
tion contribution to N given by A^fluct a- = 
-(l//?)E,{5[dctF-i(g)]/5Mj/detF-i(g). We de- 
flnc B± = mifxi^ ± mifii to establish general constraints 
on the magnitude |Ao| of the order parameter for 
s-wave pairing in the presence of population imbalance 
{Ni ^ Ni). Population imbalance is achieved when 
either -Ek,t or -Ekj is negative in some regions of momen- 
tum space. Depending on the number of zeros of -Ek.T 
and -Ek.i (zero energy surfaces in momentum space), 
there are two topologically distinct cases: (I) E^cr has 
no zeros and -Ek.-o- has only one, and (II) E'k.o- has no 
zeros and -Ek,-(T has two zeros. The zeros of -Ek.u occur 
at real momenta fc| — _B+ ± (B^ — 4mfmx|Aop)^/^ 
provided that |Aop < |B_|V(4to|TO|) for 5+ > 



and |Aop < -fi^fii for B+ < 0. The P = limit 
corresponds to case (HI), where i^k.o- has no zeros and 
is always positive. We illustrate these cases in Fig. |5| 
for iV| > N^. Notice that, the Fermi sea of the lower 
quasiparticle band is a sphere of radius k^ in case 
(I), and is a spherical shell k- < k < k+ in case (II). 
The transition from case (II) to case (I) occurs when 
A;_ — > 0, indicating a change in topology in the lowest 
quasiparticle band, similar to the Lifshitz transition 
in ordinary metals and non-swave superfluids |28i l29l |. 
The topological transition here is unique, because it 
involves an s-wave superfluid, and could be potentially 
observed for the first time through the measurement of 
the momentum distribution. 




FIG. 2: Schematic plots of i?k,T (dotted lines) and -Eka (solid 
lines) versus k a) for case (I), and b) for case (II). 

The r = momentum distributions for cases (I) and 
(II) can be obtained from Eq.^J For momentum space re- 
gions where -Ek.o- > and -Ek.-o- > 0, the corresponding 
momentum distributions are equal nk,(j = n^.-a- How- 
ever, when i?k,(T > and i?k,-cr < 0, then nk,(T = and 
rik,-CT = 1. Although this topological transition is quan- 
tum (T = 0) in nature, signatures of the transition should 
still be observed at finite temperatures within the quan- 
tum critical region, where the momentum distributions 
are smeared out due to thermal effects. Although the 
primary signature of this topological transition is seen in 
the momentum distribution, the isentropic ks or isother- 
mal kt compressibilities and the speed of sound Cs would 
have a cusp at the topological transition line similar to 
that encountered in |Ao| (see Fig.EJ as a function of the 
scattering parameter l/(kp^-^-ap). The cusp (discontin- 
uous change in slope) in ks, kt or Cg gets larger with 
increasing population imbalance. 

Next, we solve Eqs (jSJ and ^ to analyze the phase 
diagram at T = as a function of scattering parame- 
ter l/{kp^^aF) and population imbalance P = N-/N^, 
where A^± = (N^ ± Ni)/2 and k^p^ = (fc|,_^ ± fc|,|)/2. 
We perform calculations for equal (m^ = mi) and un- 
equal (to| = 0.15?7i|) masses cases, corresponding to one- 
species (^Li or ""^K only), and two-species (^Li and ""^K 
mixture) experiments, respectively. 

For T « 0, Nfiuct,a is small compared to Nq^ct for 
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FIG. 3: Plots of |Ao|, IJ.+ and /i- (units of eF,+) for m^ — 
0.15m| a) as a function of l/(fcF,+aF) when P — 0.5 and b) 
as a function of P when l/{kF,+aF) = 0. 



all l/{kF,+aF) leading to N„ w iVo,„ |30|. We define 
ep^ = kp^/{2m±) and in Fig.|21 we plot self-consistent 
solutions of |Ao|, /i+ and /i_ (in units of 6^,+) at T = 
for two cases: a) as a function of l/{kp^^ap) when 
P = 0.5 (or iV| = 37V|) and b) as a function of P when 
l/(/cF,+a_F) = (or on resonance). 

In Fig. |3^, the BCS fi± ~ (e^,! i £F,i)/2 changes 
continuously to the BEC |/i±| -^ |e6|/2, where eb — 
— l/(TO+a|n) is the binding energy which follows from 
1/.9 = Ek|rk|V(2ek,+ - efc). Since P > and the T 
fermions are in excess, all | fermions pair to form N^ 
bosons and the remaining "f fermions are unpaired. The 
amplitude |Ao| evolves continuously from the BCS to the 
BEC limit with a cusp aromid l/{kp^^ap) « —0.27. This 
cusp in |Ao| is more pronounced for higher |P| and sig- 
nals a quantum phase transition from case (I) to case 
(II), which can be detected experimentally 8]. 

In Fig.lJt), we show (on resonance) that |Ao| = (nor- 
mal phase) for P < —0.61, where /i+ w 0.64eF,+ and 
/j,_ « 0.28eF,-i-- We notice that the evolution of |Ao|, /i+ 
and /i_ as a function of P is non- analytic when \P\ — > 0, 
and signals a quantum phase transition from case (II) 
with P > to case (III) with P = to case (I) with 
P < 0. We obtain similar results when 77i| = m-i, where 
the plot is symmetric around P = 0. Therefore, this 
quantum phase transition may be studied in current ex- 
periments involving only one-species of fermions 1, 2J. 

Next, we discuss the stability of uniform superfluid- 
ity using two criteria. The first criterion requires that 
the curvature d^ilo/dAQ of the saddle point free energy 
ilo — So/f3 with respect to Aq to be positive. When 
d^ilo/dAQ is negative, the uniform saddle point solution 
does not correspond to a minimum of JIqi a-nd a non- 
uniform superfluid phase is favored. The second criterion 
requires the eigenvalues of the superfiuid density 



Py(r) = (m|iV|+mxiVi)(5,_ 



/? 



y^fcifcj-yk, 



(5) 



When at least one of the eigenvalues of pij (T) is nega- 
tive, a spontaneously generated gradient of the phase of 
the order parameter appears, leading to a non-uniform 
superfiuid phase. Notice that, Pij{T) reduces to a scalar 
PoiT) for s-wave systems. 

The uniform superfluid (U) phase is characterized by 
Po(0) > and d'^flo/dA'^ > 0, and the normal (N) phase 
is characterized by Aq = 0. The non- uniform super- 
fiuid (NU) phase is characterized by jOo(O) < and/or 
a2f7o/aA§ < 0. The NU phase should be of the LOFF- 
type having one wavevector modulation for the center of 
mass momentum of a Cooper pair near the BCS limit, 
where pairing occurs within a very narrow region around 
the lowest Fermi energy of either the lighter or heavier 
atom. Towards unitarity and beyond, we expect the NU 
phase to be substantially different from LOFF phases 
having spatial modulation that would encompass several 
wavevectors. This expectation is based on the idea that 
when the attraction between fermions gets stronger, the 
Fermi mixture becomes more non-degenerate, and there 
is a wider region in energy within which pairing can oc- 
cur, leading to a range of possible wavevectors for the 
center of mass momentum of Cooper pairs. 

As shown in Fig. ^, when m-^ = toj^, the phase dia- 
gram is symmetric around P = 0. A continuous quantum 
phase transition occurs from the NU to the N phase be- 
yond a critical population imbalance on the BCS side. In 
addition, a discontinous transition from the NU to the U 
phase of topological type (I) also occurs. 

In contrast, as shown in Fig. ^^i when m| = O.IStoj, 
the phase diagram is asymmetric around P = 0. A con- 
tinuous quantum phase transition occurs from the NU to 
the N phase beyond a critical population imbalance on 
the BCS side. Furthermore, it is found that the U phase 
has a larger stability region when light fermions are in 
excess, and that a discontinuous transition from the NU 
to the U phase occurs. The U phase also exists for P < 
when l/(fcF^+a_F) > 4.8 (not shown). Lastly, one of the 
topological quantum phase transitions (dotted lines) is 
very close to the NU/U boundary for P > in contrast 
to the equal mass case. This line indicates a change in 
quasiparticle Fermi surface topology from type (I) to type 
(II), and may lie in the U region when m^ /mi < 0.15. 

Next, we discuss superfiuidity near T « Tc, where 
Aq = 0, and derive the TDGL equation |31|. We use 
the small q and ivi —^co + iS expansion of 



where L^^{q) ~ ¥^^{q), to obtain the TDGL equation 



d 



a + h\^{x)\^~Y.-T'^NJ-^d^^ 



dt 



A(x) = (6) 



to be positive. Here, 5ij is the Kronecker delta, and in the real space x = (x, t) representation. The coeffi- 
3^k,± = (J^k.T ± 3^k.i)/2, with J^k.a = sech2(/3Sk,./2). cients are given by a = I/5 - Ek^k.+ |r(k)|V(2a,+), 



which leads to the saddle equation when 
a = (Thouless condition), and Cy = 



/3(fc,fc,C_/(TO_ek,+) 



5,,C+/2)(8Ck.+) 



_)}|rk|^ where 



C± = (Yk,T/"^T±^ka/"^i)/2, ^k,. = tanh(/3Ck,./2) and 
^k,(T = sech^(/3^k,cr/2). Notice that Cij reduces to a scalar 
c in the s-wave case. The coefficient of the nonlinear 
term is b = Ek[^k,+/(4e^,+) - /3yk,+/(8e^,+)]|^k|^ 
while d has real and imaginary parts given by d = 
lim„_,o Ek ^k,+ [l/(8el,+) + i7r<5(2a,+ - «;)/(2H]|rkP, 
where 8{x) is the Delta function. Notice that the 
damping term (imaginary part of d) vanishes for /^+ < 0, 
indicating an undamped dynamics for A(a;). 

Since a uniform superfluid phase is more stable in the 
BEC side, we calculate analytically all coefficients in 
the BEC limit where |/x±| -- |efc|/2 > T^. We obtain 
a = fli + 02 = —Vm\(2[i^ — (.\,)ap j {%t{) + Vmj^n^a^p^ 



fli + 02 
: 6i + 6; 



Vm\a\l{\%-K) 



Vm^j^y 



4 



c — y77i^aF/[87r(m| + m{)\, and d = yTO^ai?/(87r). 
Here, e labels the excess type of fermions and n^ is 
the density of unpaired fermions. Through the rescaling 
^{x) = vdA(a;), we obtain the equation of motion for a 
dilute mixture of weakly interacting bosons and fermions 



V2*(a;) .d^{x) 



2mB 



dt 



(7) 



with bosonic chemical potential /i^ — ~ai/d — 2/i+ — e;,, 
mass TTiB = d/c — m-^ + toj, and repulsive boson- 
boson Ubb — hiV/(P — ATrap/in+ and boson- fermion 
Ubf = aiV/{dne) = STraF/m^ interactions. This proce- 
dure also yields the spatial density of unpaired fermions 
given by ne{x) = [02/0? -f &2|*(a;)P/d^]/C/BF = rie - 
87raF{dne/dfie)\'^{x)\'^ /m+. Since drie/dfie > the un- 
paired fermions avoid regions where the boson field 1^(2:) | 
is large. Thus, the bosons condense at the center and the 
unpaired fermions tend to be at the edges in a trap. No- 
tice that, Eq. lO reduces to the Gross-Pitaevskii equation 
for equal masses with P = [31|, and to the equation of 
motion for equal masses with P 7^ |21j . 

Furthermore, we obtain the boson-boson gbb = [1 + 
m|/(2m|) -I- mi/{2m^)]aF and boson-fermion asF — 
AmB'me/[m^{mB + m.e)]aF scattering lengths, which re- 
duce to a^B = 2ap and ubf = S>ap/3 for equal 
masses 21, 31]. For a mixture of ^Li and "^^K, ubb ~ 
AAlap, and gbf ~ 2.03aF when ^Li is in excess, 
and aBF ~ 8.20ai^ when ^°K is in excess. For a 
better estimate, higher order scattering processes are 
needed [S^l- Since the effective boson-fermion system 
is weakly interacting, the BEC temperature is Tc = 
7r[nB/C(3/2)]^/'^/mB, where ^(a;) is the Zeta function and 
ub = {n- ne)/2. 

In summary, we analyzed the phase diagram of uniform 
superfluidity for two-species fermion mixtures (e.g. ®Li 



and '"^K) from the BCS to the BEC limit as a function 
of scattering parameter and population imbalance. We 
found that the zero temperature phase diagram of pop- 
ulation imbalance versus scattering parameter is asym- 
metric for unequal masses, having a larger stability re- 
gion for uniform superfluidity when the lighter fermions 
are in excess. This result is in sharp contrast with the 
symmetric phase diagram for equal masses. In addition, 
we found topological quantum phase transitions associ- 
ated with the disappearance or appearance of momentum 
space regions of zero quasiparticle energies. Near the crit- 
ical temperature, we derived the Ginzburg-Landau equa- 
tion, and showed that it describes a dilute mixture of 
bosons (tightly bound fermions) and excess (unpaired) 
fermions in the BEC limit. 
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